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PREFACE 


Large and small calibre guns have long been important tools of 
warfare and still are. A detailed understanding of the gas dynamics 
in a discharged gun could be of great help in the design of new guns, 
shells or charges. Such knowledge is also important to efforts to 
Suppress or detect gun smoke, gun flash, and gun blast. The work re- 
ported on in thís memorandum was prompted by a need for quantitative 
data on temperatures and densities in the gun exhaust pertinent to de- 
tection by radar or infrared scanning devices. Dr. William G. McMillan 
first directed our attention to this problem, and we are further in- 
debted to him for suggesting sources for experimental data. 

The methods are developed here for interior ballistics.  Calcula- 
tions of this sort are capable of modeling gun performance well enough 
to be an important supplement to empirical gun tests and charge firing 
trials. Іп addition, some understanding of gun flash, blast, and smoke 
have been gained by extensions of these methods. 

We are greatly indebted to the personnel of the U.S. Naval Weapons 
Laboratory, Dahlgren, Virginia (principally Dave Bowen, with the coop- 
eration of Frank Kasdorf) for the collection and interpretation of gun 
characteristics and firing sequence data used for inputs to and compari- 
son with our calculations. 

Helpful discussions with Alan Shapiro here at Rand are also ac- 
knowledged as is the prior interest in this subject by Jack Craig. Ме 
would also like to thank Wakichi Asano for his valuable programming 
assistance. 

This work represents a progress report, since the authors can 
foresee further investigations of charge characteristics, energy losses 
in the barrel, and other non-ideal effects, as well as a more detailed 
hydrodynamical model of the smoke and flash. Expressions of interest 
in such extensions and criticisms of the work presented here are wel- 


comed. 
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SUMMARY 


This paper considers the gas dynamics of various large naval 
guns, including the interior ballistics (charge burning, propellant 
gas expansion and shell acceleration) and the nature of the. gas expan- 
sion and air shock beyond the muzzle. Solutions are accomplished by 
a numerical program which integrates the partial differential eguations 
appropriate to the gas dynamics in one space dimension, Тһе eguations 
of state for the gas and the air are treated in some detail, as is the 
nature of the charge burn, the barrel friction, and the shell inertia. 
The dynamics of the partially burned charge gases which sometimes 
flash burn after ejection from the muzzle are treated within the limits 
of expansion into various conical geometries. Comparisons are made 
with observations of both the interior ballistics and of the smoke, 


flash and blast from the exhaust. 
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I. INTRODUCTION 


This paper introduces an accurate numerical method for obtaining 
the hydrodynamic properties associated with the propellant gas as it 
moves down a gun barrel and exhausts to the air. Previous analytical 
solutions to this problem have found it necessary to make a number of 
simplifying assumptions and approximations, some of which limit the 
accuracy and range of applicability of the results. Even the very 


14) takes no account 


recent treatment of high speed guns by A. E. Seiga 
of the nature of the propellant and simply assumes that it burns in- 
stantly into an ideal gas with constant specific heat. Further, Seigal 
assumes that the shell moves without any frictional or heat losses. 

By far the most extensive treatment of interior ballistics is by 


(2) 


J. Corner. He summarizes all important work in the field before 
1950. Тһе "isothermal" solution is the simplest and it assumes that 
the propellant burns according to a specified law into a gas which ex- 
pands isothermally during burning with a constant specific heat ratio 
which is increased from that of an ideal gas to account for heat losses. 
The shell mass is increased to account for frictional resistance. 

Generally, in Corner's theory, in order to arrive at the correct 
peak breech pressure, muzzle velocity, pressure-time history, and shell 
ejection time for a given gun, the solution is made to depend on three 
propellant parameters, the burning law, the total energy released, and 
the rate of burning. Corner describes a more detailed solution which 
takes account of the temperature dependence of the gas and includes а 
number of other small corrections. It leads to very complicated ana- 
lytical expressions which depend basically on a propellant burning law 
parameter and a burning rate parameter. Corner shows a few cases where 
good agreement with experimental values can be obtained by proper ad- 
justment of the parameters. Each gun must be treated separately and 
the solution is very sensitive to the various adjusted parameters. 

The treatment we present here takes good account of many of the 
factors in the problem which previously suffered from simplifying 
assumptions. Using a full partial differential formulation of the gas 


dynamics simplifies some factors which have only a minor effect on the 


results and reduces interior ballistic problems to a single free param- 
eter, namely the propellant burning rate. This parameter is obtained 
for one gun by fitting to experimental data and the solution for peak 
chamber pressure, muzzle velocity, pressure-time history, and shell 
ejection time is then shown to be accurate to the order of ten percent 
for several different guns. 

Previous descriptions of gun flash and smoke have almost entirely 
been qualitative in nature and no rigorous mathematical treatment is 
known to exist. The hydrodynamical solution which we present, while 
limited in its thoroughness by a number of simplifying assumptions, 
does reproduce quantitatively a number of properties observed in actual 
guns, such as, the pressure as a function of time and distance from 


the muzzle and the shock reheating of the gas outside of the gun. 


ІІ. INTERIOR BALLISTICS 


In the results presented here are examples for six different 
naval guns ranging in size from 5"/38 caliber to 16"/50 caliber. The 
physical characteristics of these guns are listed in Table 1.) The 
table includes such parameters as chamber volume, propellant weight, 
shell weight, and shell travel distance. All the propellants used in 
these guns are essentially identical. They consist of roughly 90% 
nitrocellulose, 2-3% butyl stearate, 1% ethyl centralite, 1% lead 
carbonate, 1% potassium sulphate, 3% volatiles, and 1-2% moisture. ` 
They come in the form of cylindrical grain roughly of diameter 4, 
and length % = 24, with roughly six web perforations of diameter 
wel Values for these quantities are also given in Table 1. 

In this analysis of interior ballistics, we have made the follow- 
ing simplifying assumptions: (1) each gun is unworn so that no gas 
leaks between the shell and the barrel; (2) the only losses in the 
system are frictional in nature; (3) the pressure drop needed to over- 
come the "skin friction" of the barrel in moving the gas down the 
barrel is negligible; (4) gun recoil and the energy expended in shell 
rotation are negligible; (5) the chamber can be replaced by a cylin- 
drical continuation of the barrel with the same total volume; and (6) 
the propellant grain distributes itself uniformly between the breech 
and shell during its burning. These assumptions can be shown to be 


a CO 


consistent with the general accuracy of our metho Other conditions 


necessary in formulating the solution are discussed below. 


MATHEMATICAL FORMULATION 





The eguations which describe the interior ballistics part of the 
calculation are one-dimensional since the cross section is assumed to 
be uniform along the length of the gun. Тһе eguation of motion іп 


Lagrangian form is 


82 --2 PQ (1) 
_ St) 


where u = u(x,t) = Sc: is the velocity of the gas in meters/milli- 
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second; X(x,t) is the Eulerian coordinate which gives the position, 
at tíme Е, of a gas element that was initially at position x, the 

| х 
Lagrangian coordinate; m = m(x,t) = 2 p(x,t)dx is the mass of the 


d 2 : Р P 3 
gas in megagrams/meter ; р = p(x,t) is the density in megagrams/meter'; 


3 3 6 
P = P(x,t) is the pressure in jerks/meter” (1 jerk 101 ergs or 


1 ik = 10“ bars), and 


2 
Q = Q(x,t) = > | = р рай | for 2720 (2) 
? V at V at ot 
oV 
and = 0 for SE > 0, 


(5.6) 


is an artificial viscosity pressure introduced to remove mathe- 
matical discontinuities which would otherwise occur whenever compres- 
sions lead to shock waves in the gas. Typical values appropriate to 
this problem are с. - 6.0 апа с, = 0,5. Ап is а mass element апа 

V = V(x,t) = l/p(x,t) is the specific volume of the gas in шанага) 
megagram. 


The equation of energy conservation is 


ӘП =- (P+ +3 (Р) (3) 


where E = E(x,t) is the internal energy density of the gas іп jerks/ 
megagram and B(P,f) is a source term which represents the rate at 
which energy is being added to the gas by the burning of the propellant. 
B can be written as a function of pressure (P) and of the fraction of 
propellant already burned (£), 

For most propellants the burn rate in energy per unit time (B,) 
is known to be directly proportional to the exposed propellant surface 
area (S) and the gas pressure on the propellant (P), or 


B, = cSP (4) 


where с is a constant. This relationship has been verified for nitro- 


(8) 


cellulose, although at pressures near atmospheric P should be 


replaced by a + ЫР, where a and b are constants, 


For the cylindrical grain described earlier, neglecting the burn- 
ing of the ends and the webbing as a first approximation, we can write 
5 as a function of Ғ and other known constants. For the fraction of 


total propellant burned as a function of time we can write 


2 2 
f= (0) -1- EI 1. 210 (5) 
ndo 4/4 do 


Now we assume that the propellant of total weight M EBD consists of n 


grains each of weight UN and density D 50 


= = ша 2 
pes пи = ne Vi зар, Tdo 0/4 (6) 
and 
4 Y x d 
S = ns, = midl = —ЁЕОР. (7) 
рдо 
so, using Eg. (5), 
aW ro 3 aW r 1/2 
Bo ote E 52 c Prop (1-8) p (8) 
Е р do 
Pado n 


Thus we can write B(P,f) in energy per unit mass per unit time as 


tan ee | (9) 
AW гор 1 4 
where s = c SEE = .01 m /Mg-ms/d,. Тһе numerical value .01 has 


n 
been estimated, independent of Moon and NE from burning rate data (9) 


and then modified to give a better fit to the experimental data, as will 
be discussed later. Also, we assume our nitrocellulose propellant has 
a calorific value, E rop & 1000 cal/gm - 4.2 jerka/Me: ©) 

The equation of motion of the shell is obtained by using the ve- 
locity and pressure of the gas at the shell boundary as it moves down 


the barrel: 
du(x >t) 
Mee oa = P(x pt) A- Fer (10) 


where A is the cross sectional area of the barrel and Fr, assumes that 


-7- 


all the energy lost in the motion of the shell is caused by а friction- 
al force acting between the shell and the barrel which is directly 
proportional to the pressure on the shell (appropriate Жос a fluid on 
a deforming plastic seal or metal shell): 


Fer = fr • P(x pë) A. (11) 


Eq. (10) can be rewritten as follows: 


du РА КОЕ: УЕЛ På (12) 
sh eff 
M 


_ sh Е Ек 
where M E re 





Тһе fraction of kinetic energy of the shell going into frictional 
energy loss is thus fr/(l-fr). Corner suggests that this fraction is 
roughly twenty percent. (10) 

Other models for the dissipative force were tried but seemed to 
ba much poorer in terms of physical soundness and resultant solutions, 
One such alternate scheme assumed the frictional force, Еа. (11), was 
a constant and another subtracted energy from Eg. (3) as a function of 
the kinetic energy of the shell, 

In an effort to derive an appropriate caloric eguation of state 
for nitrocellulose we have relied on the similarities obvious in the 


u we were unable 


Chapman-Jouguet adiabatics for other explosives, 
to find any experimental detonation data for nitrocellulose. In Fig- 

ure 1 we have plotted Y- 1 vs р for TNT with loading density ро = 

1.50 d Sé and for Pentolite with pg = 1.65 sm/em>. (13) By fitting 
а smooth curve to this data in the region of interest, namely, р < 


„5 sem we obtain 


2 
Y(p)-1 ~ 209 + 202 (13) 


.3 sb ae 


Nitrocellulose propellant generally has a loading density Po 5.7 


3 
gm/cm . This lower initial density prompted a somewhat arbitrary 
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Fig. 1—Caloric equations of state 
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modification to the fit: 


„09 + gå 


„3+ p 


(14) 


using the fact that the value of Y(p)-l in the high density limit 
(р>>1 ale) has been shown to be linearly related to loading density. 


(14) 


The formulation of this result by Miller, based оп explosives, shows 


that we can write 


У(®)-1 = + бро/р, (15) 


where ро is the loading density of the ехріовіуе, Po is its crystalline 
density, and G is a constant determined by the Morse atomic parameters 
associated with each of the elements composing the explosive, An 


independent determination by Fickett and oba, с») 


using various thermo- 
dynamic relationships, along with experimental data on the thermal 
equation of state and detonation velocity for various explosives, leads 
to empirical results for У(о)-1 similar to Еа. (15), The transition 
from the high density region to the low density region of the gas is 
obtained by using Еа. (14). 


Thus our caloric eguation of state is 


2 
P = Ep[y(p)-1] = Ep EEF (16) 
+ p 


For low densities, р<<,1 čen Eg. (16) describes an ideal gas of 
У- 1.30, The effect on the results of varying the density dependence 
of y(p)-1 has been investigated and is reviewed in the Results and 
Discussion Section. 

The thermodynamic description of the gas is completed by using a 


thermal equation of state of the form 


T = (WR) (Р/р) / (1+ Bp) (17) 


where R is the Rydberg gas constant, M is the average molecular weight 


-10- 


of the propellant gas (about 30 gm/mole), and B(p) is a "covolume" 
corraction (16, UD the ideal gas law, which is roughly linear in 
density and about 1,-2. at р = .5 sue. depending on the exact form 


of the correction used. 
NUMERICAL SOLUTION 


A solution for the hydrodynamic functions in this problem results 
from numerical integration of the two equations of motion, Eqs. (1) 
and (12), the energy equation, Eq. (3), and the caloric equation of 
state, Eg. (16). These can be written in difference equation form as 


follows: 


m+1/2 _ n-1/2 АЕ роп n-1/2 ол-,24 (18) 
Se вэ. mm, ІР М2 7 74. HO 8%ы/2.7 41/24 
„1/2 | т-1/2 e A př (1-fr) T 

j max ~ "j max Мо å max 22 (19) 

n+l cH leg n+1/2 
£l Eus Нені * Tu 5 4 Wy aye 7 Yp An 

e 1/9 (20) 
+s Ру) (1- [E улі, 

nli _ ntl n+1 nri n+1 
Pea 7 Ej-1/2 P3-1/2 C409 + (р; 15 å (0. 1/2 A: xen 

with 
27 = x o { EV? (22) 
n+1 үл 2 
qur? _ "mI T (Vi-1/2 7 V4-3/2 
j-1/2 ntl n+l/2 2 
v 1/2 * Yi- 1/2) (At у 
yt y? 
p Za Pra и 
n+1/2 (23) 


ntl 
3-1/2 * jap? ^E 
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n+l n Ё ntl j 1-1 
for Vi-1/2 < Уі 1/2: with Vi 179 = Тэн (24) 


where n 2 land 1 Í 1 Í j шах. Тһе mass elements Am, үүр are fixed 


initially and remain unchanged throughout the solution. Notice that 
in writing the source term of Eg. (20) we have used the pressure and 
energy density of the breech. The pressure and energy density are not 
actually uniform throughout the gas, but as discussed later, these 
approximations have a negligible effect on the results. | 

Unit increments in the superscript n represent step-wise incre- 
ments in time and unit increments in the subscript j represent finite 
increments in the mass, corresponding to a spatial zone, where j = I 
is the zone at the beginning of the chamber and j = Tang is the zone 
adjacent to the shell. The integration is performed by advancing in 
time (n) with the step size determined by stability conditions which 
depend on the local sound speed in the gas and the artificial viscos- 


(18) 


ity. Then for each step in time, we obtain а complete hydrodyna- 
mical description of the propellant gas by advancing in mass zone from 
j=1ltoj= Тағы These eguations yield а self-consistent solution 
for u, E, p, P as functions of time and position, with T determined 
independently from p and P. 

By using the loading density and an initial energy density of 
the propellant assuming atmospheric temperature, we specify gas condi- 
tions in the chamber. Тһе constants fr, which determines the effect- 
ive mass of the shell, Е тор? the energy/mass released from the pro- 
pellant, and s, the rate of propellant burning must also be specified. 
From previous discussion, we know fr = .2, E rop = 4.2 jerks/Mg, апа 
s = .01/d,. These values are used initially and then modified to ob- 
tain the best possible results. Тһе effect of varying these parameters 


is discussed in the results. 


RESULTS AND DISCUSSION 


In Table 2 we present for each of the six guns which we have 


studied the best fit values for the three guantities, Maus the shell 
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muzzle velocity, Pax’ the peak breech pressure and iu the shell 
ejection time, along with the corresponding experimental values. In 
Figure 2, we present calculated and experimental pressure-time curves 
for the 5"/54 gun. А11 solutions were obtained using fr = 0, and 

E rop = 4.2. We used s = 1.48 for the 5"/54 gun and then computed 

si = „01/47 as an estimate for the other guns. The actual value of 

s used along with s“ are shown in Table 2. Тһе results show agree- 
ment between the best fit values and experimental values to be better 
than 5Ž in most cases and better than 10% in all cases. 

In Table 3 we present a summary of the efficiency of guns. There 
are two standard measures of efficiency. The first is called "piezo- 
metric efficiency" and is defined as the average pressure necessary to 
accelerate the shell to the observed muzzle velocity divided by the 


peak breech pressure. We compute the average pressure P as follows 


from Eg. (13), with fr = 0, 


42 
E T 2 (muz ) 
PA= Mon а = M, TR | (25) 
tot 
2 
_ 1/2 Mu 
or P= Ew AES (26) 
tot 
where Re GE is the combined chamber length, Len? and shell travel length, 


5, 09) 


. Piezometric efficiencies are usually about .4 - We see 


Leray 
from Table 3 that the efficiencies for these guns are somewhat lower 


than this. 


' and is defined as 


The second measure is "ballistic efficiency,' 
the ratio of the kinetic energy of the shell at muzzle velocity to 
the energy of the propellant. We shall calculate this efficiency in 
two ways: first we shall use the full energy/mass of the propellant, 
= 4,2 jerks/Mg, and second we shall use the energy/mass actually 
"prop (т.ј 


values along with B run (Toi) > are shown in Table 3. The typical bal- 
listic efficiency of 1/509] lies roughly between our two sets of cal- 


E 
prop 


burned at the time of shot ejection, ). Both these efficiency 


culated values. The maximum possible ballistic efficiency for these 


guns, obtained if the charge is completely burned before the shot moves 


appreciably, can be shown to be roughly 1/2, 09) 


Breech pressure, P (kilobars) 
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Fig.2—Pressure-time history for 5"/54 gun 
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Іп order to interpret the significance of these results, we must 
discuss the changes caused by varying the parameters in the problem. 
These changes were obtained by keeping fixed all parameters but the one 
whose affect we are studying, and then obtaining solutions for several 
values of this parameter. Changing fr from .0 to .25 decreases Üiz 
by less than 4% and increases Eins by about 25%, with the changes in 
Up and LER roughly linear with the changes of fr. Changing Е тор 
from 4.2 to 3.0 jerks/Mg without changing в, decreases Ша by about 
3% and increases Ры about 1%. One reason that this change has so 
little effect is that for = = 4.2 jerks/Mg the propellant is only 


rop 


about 60-75% burned at T ., and lowering the value of E to 3.0 
ej ргор 


jerks/Mg simply increases the fraction of propellant burned at Taj“ 


Lowering E rop still further would then begin lowering Be апа as 
quite strongly. The fact that the nitrocellulose propellant we have 
used is only partially burned at Tos indicates that it most likely has 
a calorific value less than 1000 cal/gm. Modifying the form of y(p)-1 
from Eq. (14) so that ік goes from the minimum of .3 at low density to 
a maximum which is 1.2, increases U iuz by about 10% and s by about 
20%. The shape of pressure-time history is somewhat less sensitive to 
changes іп y(p)-1 but if the maximum value is raised significantly 
above 1.0, the peak in the curve occurs sooner than 1/3 Sg and if the 
maximum value is much lower than 1.0, then the peak occurs at roughly 
1/2 тар The results аге fairly dependent оп ү(р)-1 but we have chosen 
an expression which we believe to be most accurate. 

After we have established values for fr, Е гор? апа ү(о)-1, as dis- 
cussed earlier, the guantity s is our only variable parameter. We study 
the effect its variation has on ЦЭЭР and шин for each of the guns. 
Using the fixed conditions fr = 0. and У ар = 4,2, with у(р) given by 
Eg. (14) and with B(P,f) given by Eg. (9), a curve of Uz Versus 8 is 
shown in Figure 3 and a curve of Базы versus в is shown in Figure 4. 

Combining Figure 3 and Figure 4 and plotting peak breech pressure 
versus muzzle velocity in Figure 5 it is seen that as the burn rate 
increases the pressure rises much more rapidly than the velocity. This 
reflects the fact that as the burn rate increases the pressure builds 
up faster to a larger peak value, while the amount of energy put into 


the gas approaches the total energy in the charge. 
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Fig.3—Muzzle velocity versus burning rate 
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We have also investigated using burning laws of the form B(P,f) = 
s P^(1-£Y) where .8 < a < .1 and 2 < y 5 10 and have found that Мач 
and Es are not particularly sensitive to о and y as long ав в is ad- 
justed to keep the burning rate roughly the same.  Às mentioned earlier, 
we have also made calculations using individual zone pressure and energy 
instead of breech pressure and energy, and these yield the same results 
as long as s is adjusted slightly to keep the overall burning rate the 
same. In addition we have investigated burning the propellant accord- 
ing to Eq. (10) only in the chamber part of the gun. We have found that 
for a fixed value of s,burning only in the chamber vs. burning all the 
way between the breech and the shell leaves Bine approximately the same 
and decreases i un by about 102. 

Às a further study on the nature of the burning, we have initially 
ignited oniy the first 20Z of the propellant. This is done by setting 
the initial energy density of all but the first 20% equal to zero and 
then letting the hydrodynamical calculations propagate the ignited gas 
down through the remaining propellant. The result of this is to create 
shock waves in the gas which propagate between the shell and the breech. 
The magnitude and time dependence of the resultant oscillations in breech 
pressure depends critically on the fraction of propellant initially ig- 
nited and also on the burning law for the remaining propellant. We be- 
lieve that this initial burning process thus accounts for the oscilla- 
tions which are observed in the actual pressure-time histories. However, 
since we do not know the precise experimental details of the ignition 
process, we have retained the original assumption of uniform burning as 
a good first-order approximation, consistent with accuracy of our other 
assumptions. 

Às a final study we have investigated the effects of mass ablation 
and turbulant heat transfer on our results. These energy loss effects 
are treated independently of the frictional losses and require modifi- 
cation of the equations of motion and energy. In the Appendix we give 
a detailed description of the new equations and show that the results 


are only slightly affected by the changes. 
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ІІІ. FLASH АМ) 5МОКЕ 


We now treat the propellant gas as it exits from the muzzle. Тһе 
gas composition is such that it usually forms a combustible mixture 
with the surrounding air and а momentary flash of bright light results 
from a shock-initiated burn. Based on our survey of available experi- 


(3) (20) 


mental data and on previous studies, the flash phenomenon can 
be divided into basically three regions. The first two regions are 
muzzle glow, a small region of low luminosity at the muzzle, and the 
primary flash, a high-luminosity region which is separated from the 
muzzle glow by a dark region. These regions are usually discernible 
only for small guns. For the large guns which we are considering, the 
third region, or secondary flash, is a highly luminous continuation of 
the primary flash and occupies most of the gas volume. 

The cause of the flash phenomenon has been explained gualitatively. 
Basically, the gas expands and cools very rapidly upon leaving the 
muzzle, and conseguently the luminosity diminishes and the dark region 
results. This gas is then recompressed through a decelerating shock, 
which raises the temperature back up to roughly the muzzle temperature 
and forms the luminous primary flash which is relatively small and 
usually lasts for a few milliseconds. Тп the meantime air has mixed 
with the unburned hydrogen and carbon monoxide in the muzzle gas, form- 
ing a combustible mixture. The exact composition of the muzzle gas is 
highly dependent on such parameters as the initial loading density of 


(12) 


the propellant and the final density of the gas, -besides being 
dependent on the original propellant composition. ТЕ the temperature 

of the recompressed gas is above the ignition temperature of the com- 
bustible mixture or if a piece of burning propellant penetrates the gas 
cloud, then it will burn and form the large secondary flash. Тһе occur- 
rence of flash is often inhibited by the presence of certain substances, 
such as potassium sulphate, which act as suppressants. If the unburned 
propellant gas is only partially ignited or not ignited at all, then 

it appears as smoke outside of the gun. 


In order to form a hydrodynamical model of the flash phenomenon 


we will make a number of assumptions which simplify the calculations: 
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(1) the muzzle дав emerges into а radially symmetric conical section, 
which is initially a full hemisphere and then steadily decreases to a 
very narrow angle cone; (2) the emerging muzzle gas displaces the 
surrounding air and no explicit account is taken of any mixing between 
the two; (3) the unburned propellant remains uniformly mixed throughout 
the gas, but the chemical composition of the gas itself is not deter- 
mined or used in the calculations; (4) the emerging propellant gas 

obeys Eg. (14) and the surrounding air is treated as an ideal diatomic 
gas of specific heat ratio y = 1.4. Although these assumptions are 
guite restrictive, we still expect to gain certain guantitative informa- 


tion about the nature of the flash phenomenon. 
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MATHEMATICAL FORMULATTON 


The numerical solution now separates into two parts. Тһе interior 
ballistics are still governed by the previously described eguations of 
motion and caloric eguation of state. For the expansion into air we 


have the equation of motion 


ou 29. 


ar ‚Sm, rn) en 


ot 


where r, the radial distance from the muzzle, now replaces X, the 


linear distance used іп Eg. (1), M = le p(x st) x" аг = ын Р үз із 
the total mass іп an Fth of а sphere where F < 1, De = Мат = Ep 
270 


r is corresponding mass per steradian. Ме can write F = sin 2 


where 9 is the angle between the axis and the boundary defining the 


cone into which the gas is expanding, and 9% is the modified artificial 


viscosity 
2 
F px 4 Әк 2 ot Әк 
Vr Vr 
2 (28) 
where C3 and с, аге constants and Am, is a mass element, 
The energy eguation outside the gun becomes 
OE ov 
ДЕ (Р + Qg) At + SR (29) 


where SR is a source term which represents the burning of the gas. 
We shall arbitrarily assume tbat SR can be written as a function of 


gas temperature T in a form similar to Eq. (9), 


1/2, 


SR = с(1- £) (30) 


where c is a constant which will be adjusted so that the gas will be 
completely burned about 10 milliseconds after first exiting from the 
muzzle. Т is obtained by using the thermal equation of state, Eq. (17), 


which reduces to the ideal gas law at the low densities outside of 
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the gun, Тһе air surrounding the gas does not burn and so we set 


SR = 0. in that region. 
Тһе caloric equation of state for the gas is given by Eg. (16) 


and the eguation for air is 
Р = .4Ep (31) 


Тһе approximation of air as an ideal gas is valid for the relatively 


low temperatures and densities which occur in the flash region, 
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NUMERICAL SOLUTION 


We can now write Eg. (27) to (31) in difference eguation form 


2 
n n 
QUU. m-U2_ ший бй: A da n-l/2 _ y n-l/2 39) 
j j Amp j+1/2 1-1/2 Р 41/2 F. 1/2 
ml _ 1, mti n+1/2 mH n 
E12 = Eja "3053/2 + ”4-1/2) + % GER 1 01-12: 7 “141127 
n n 
+ SRY ayo At (33) 
mtl _ n+l nri 
Pi-1/2 755 Еу руа 5-1/2 (34) 


where, in the gas region, 


1/2 


n = n 
SR, 172527, 1/2 Он E B. 1/2 Etot (35 
and where 
n+l ya 
«+1/2 20 Ames 1/2 SCH 1/2 7 -1/20 
¥i-1/2 Ары п . ARIPO. A BET? 
1-1/ "INE Or Aat 2)2 АС H2 
(36) 
nti n 
+ 20, Алы) туз 1У74-1/2 7 У1-1/2 | bey ev 
Che + Yigg (еп Лу Az PL 
and 
: se 
A"gi-1/2 7 (37) 
j-1/2 
with 


өші, „ml | 2 
аа, 2014-1214 PR Ge > Do/2 (38) 
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2 
ntl 
JN Dei for ge < Do/, and SE eng 
2 
and 
пк 3 n+1, 3 
nil Ap [Cr - (r 1) 1 ntl / 
E жк Ме / E 
всс т RT pes 
ntl, 3 2 
pU A а 
3 4 ©? GRECI sat 
D 
n+l о 
ES cm 


where D, is the gun barrel diameter. 
The solution is obtained by defining the origin as the center of 
the muzzle and supplying as initial data the velocity, pressure, den- 


sity, and temperature of the propellant gas as a function of negative 


distance from the muzzle, үй s X" - R for 0. < x7 < В and 153 
j tot j tot 
= Te at the time when the shell reaches the muzzle. Тһе numerical 


integration uses Eqs. (18), (20), (21), (22), (23), and (24) for the 
plane geometry region т} > Dol» where Eqs. (38) and (39) give the 
appropriate geometrical factors for the transition between the two 
regions. 

The initial conditions for the air outside the gun are obtained 
by assuming standard temperature (T = .03 x 10“ 5 
(P = 9.5 х 10? jerks/m”). 


These calculations were carried out for the 5"/54 gun as а repre- 


K) and pressure 


sentative case and one for which we have comparable experimental даса. 
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RESULTS AND DISCUSSION 


A survey of several dozen firings o£ the 5"/54 gun, as recorded 
on 16mm film taken at 96 EE shows flash and smoke character- 
istics which аге in good agreement with the general description given 
earlier. The gas initially emerges into a full hemisphere (0 = 1/2, F 
= 1/2), after about 10 milliseconds it emerges into a cone of angle 
0=11/3 (F=1/4), and after about 20 milliseconds it emerges into а 
cone of angle 0= 1/4 (F=1/8). After about 50 milliseconds the gas 
is emerging straight forward, which it continues to do until it is 
completely exhausted about 100 milliseconds later. 

In the firings where there is no flash, the initially emitted 
gas forms a large cloud of smoke roughly two meters long and two meters 
wide which drifts away from the muzzle slowly, and through this is 
shot the gas which emerges at a smaller angle and forms a second cloud 
of smoke further downstream. In the firings where there is a flash 
it is possible to see the dark region of expanding gas just outside 
of the muzzle, and then about 1,-2. meters downstream the shock forms 
which initiates the flash burn within the first 10 milliseconds. The 
extent to which the flash spreads throughout the rest of the gas de- 
pends guite a bit on the amount of flash suppressant in the gas. Some- 
times the flash is only a small region of luminosity which extinguishes 
itself very soon, and sometimes the flash consumes all of the gas and 
is very bright, extending over several meters. 

We have obtained by hydrodynamical solutions for the gas and air 
outside of the gun for the cases F= 1/2 (9 = 1/2), F = 1/4 (8 = n/3), 
and F = 1/8 (0="/4). To show the formation of the shock, we have 
plotted gas temperature versus radius at a fixed time for each of the 
three values of F (see Figure 6). Note how the gas cools very rapidly 
as it comes out of the muzzle and then at a distance of 1.5 - 2, metars 
the temperature raises back up to its initial value as a compression 
shock is formed. The unburned gas is then ignited and flash results, 
Note that the location of the shock region does not change greatly 


with a change in F. The fact that the shock distance at a given time 


Temperature, Т (105 °K) 
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Fig. 6— Temperature-radius profiles for 5/54 gun 
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is greater for smaller F shows that the gas moves faster when confined 
to a narrower cone. In addition to temperature, our calculations give 
the gas pressure, density, velocity and internal energy as a function 
of position and time, 

In Figure 7 we have plotted the gas pressure (assuming F = 1/2) 
at a fixed distance from the muzzle (r = 2.m) versus time, Note the 
sharp rise in pressure as the gas first arrives at 2.m and the gradual 
decrease back to atmospheric pressure. This pressure-time history is 


(3) 


in good agreement with corresponding experimental data, except the 
experimental pressure curve usually drops somewhat below atmospheric 
and then forms a second peak of roughly the same magnitude as the 
first about 10 to 15 milliseconds later. Based on our experimental 
data, this second peak is due to the flash burning of the gas and it 
does not occur if there is no flash. The fact that we do not reproduce 
a second peak is probably due to the fact we do not burn the gas in 
the air properly, which would require knowledge of the exact chemical 
composition the gas and its mixing process with air. Also there are 
small instabilities in our numerical integration procedure, caused by 
the transition from plane geometry to conical geometry, which became 
larger with increasing time. These cause spurious oscillations in the 
pressure profile which would most likely obscure a second pressure peak 
even if it was caused by our burning procedure, 

In Figure 8 we have plotted peak pressure versus radius using our 
calculations for the three angles 9 = 1/2 (F = 1/2), 0 = m/3 (F = 1/4), 
and 0-т/4 (F = 1/8), along with the actual experimental data measured 


(3) 


at these angles." “Note that we have made the assumption that the 
pressure at a given angle is obtained by using the value of F corres- 
ponding to that angle in the numerical integration, Ме feel that this 
is a rough approximation to the actual pressure distribution, Although 
the calculated pressures are greater than the corresponding experimental 
values by about 50%, they all lie on straight lines with almost the 

same slope, thus correctly reproducing the adiabatic expansion of the 
gas as it moves away "from the muzzle, Тһе fact that we have had to 


rely on separate calculations using different conical angles in order 
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to simulate the angular dependence of the pressure points up а major 
limitation in our one-dimensional conical calculations. 

A final fact which comes out of our calculations is that the time 
for the pressure peak to travel from a radius r = 2m to r = Dm is 
roughly 6ms, which is in good agreement with the experimentally meas- 


ured time. 


-33- 


IV. CONCLUSIONS 


We have obtained a good quantitative model for the interior bal- 
listics of large naval guns. Our results are limited by a number of 
simplifying assumptions but they are valid within these limitations. 
Our major improvement over previous treatments is the use of a more 
realistic equation of state for the propellant gas. We have shown how 
sensitive the interior ballistics are to a number of parameters. With 
more detailed experimental data we improve these parameters somewhat 
and also investigate the replacement of the various ideal assumptions 
by the actual physical process. 

The flash and smoke are much more difficult to model in detail. 
Although our treatment accounts fairly well for the pressure distribu- 
tion and shock reignition of the gas, we need a much better knowledge 
of the mixing of the gas and air and the flash burning mechanism. A 
number of the mathematical features of the hydrodynamics need to be 
improved. The use of two-dimensional calculations is necessary to pro- 
vide for the angular dependence of the pressure. Alternate methods of 
treating the transition region between plane and conical geometry should 
be studied in order to improve the stability and accuracy of the calcu- 
lations in that region. 

We have established the validity of our calculations for large 
naval guns, but in order to show their general validity they should 


be applied to a variety of other guns. 
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APPENDIX 


A study of some gun barrel effects which have a small but finite 
effect on the gas dynamics can be made by generalizing the Lagrangian 
equations of motion to include the effects of heat and mass transfer 
due to ablation and gas turbulence processes during the interior bal- 
listics. 

The equations of motion are recast іп an "almost-Lagrangian" form 
which allows each of the zones to increase in mass as wall ablation 


(21) 


occurs. Following the derivations of Crowley, we obtain an equation 


of continuity 


Speed x 39V а а = 
eg C ae bns (A-1) 


where m is the mass flux (mass/area/time) and S, is the surface area 


per unit mass. The momentum equation becomes 


du _ ER 2 E = E 
scm on (P+ 0) +m Ss, (u, u) Tom | (А-2) 


where чу is the velocity of the mass entering а zone and pe is the 


shearing stress at the wall. The energy equation becomes 


= =- (р + 0) zZ +m S Г(а-а 72/2 +E,- E] (A-3) 


+ T8 Ju] +H 
wm 


where the second term on the right hand side gives a rate of energy 
addition associated with the mass addition and B is the specific in- 
ternal energy of the mass entering a zone. The third term comes from 
work done on a mass element by wall surface forces, and the fourth 
term, H, is the rate at which heat is supplied to an element by turbu- 
lent heat transfer. | 


The following expression for the mass flux is used: 


Си (YPu) / (ү-1) Cpu 


Xx T 4x PDE T 
E + n(PV+E) n + Е / (PV+E) 


(A-4) 


n= 
q 


where Cy is the dimensionless coefficient of heat transfer (the Stanton 
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number), E, is the sum of the enthalpy per unit mass required to heat 
the wall material to vaporization temperature, the specific heat of 
fusion and the specific heat of vaporization, and n is the turbulent 
transpiration coefficient. The source term, representing the heat lost 


from a zone due to the turbulent ablation process, is written as 


H=- а Sm =- Cp (РУ+Е)и5 (A-5) 


The shearing stress can be written as 


ы 2 _ 2 х 
Сэгс Ce pu = On pu (A-6) 
where Cc is the dimensionless coefficient of skin friction and Ca = ic, 


by the Reynolds alder 


To obtain а numerical solution we now put Eqs. (А-1) to (A-6) in 


difference equation form. Using 


n 2 n_n xod m п = 
dod =D (RG - XP) “р, Vio gai; gs (A-7) 
n n-1/2 
C. p. u. 
. Og P4-1/2 53 _ 
EDT > : п n : a) 
ME А 112 
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41/2 j-1/2 n,_172 
noe n n n n 2 
== Cy 0513/2. (55:3 /2V 424 723844 72? (А-9) 
: ae gn+1/2 ‚ 
3 Ші. 1/2 
2 
n-1/2 14,5, 1 , n-1/2 n-1/2 п-1/2 = 
and Tu, 2562 Oju 594-120 OG 0) CASTO) 


and substituting into Eqs. (A-2) and (A-3) we obtain 
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